Phase diagram of the one-dimensional extended attractive Hubbard model for large 

nearest-neighbor repulsion 
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We consider the extended Hubbard model with attractive on-site interaction U and nearest-neighbor 
repulsions V . We construct an effective Hamiltonian -ffe// for hopping t «V and arbitrary U <Q. 
Retaining the most important terms, -ffe// can be mapped onto two XXZ models, solved by the 
Bethe ansatz. The quantum phase diagram shows two Luttinger liquid phases and a region of phase 
separation between them. For density n < 0.422 and U < —4, singlet superconducting correlations 
dominate at large distances. For some parameters, the results are in qualitative agreement with 
experiments in Bai-ajK^jBiOs. 
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I. INTRODUCTION 

BaPbi_a;Bi2;03 and Bai_a;Ka;Bi03 are suBerconduG=. 
tors with transition temperatures Tc near 13KEI and 30Ktl 
respectively. In spite of the fact that these systems do not 
contain Cu and are non-magneticEl, they have important 
similarities with the cuprate superconductors, like the 
perovskite structure, relatively high and low density 
of states at the Fermi levefl. While the electron-phonon 
interaction is very important in these, materials, as shown 
by inelastic neutron measurementso, and the existence 
of displacements of O atoms and lattice distortionsElQ, 
several other experiments suggest that the pairing mech- 
anism is, at least in part, of electronic naturm. Exci- 
tons combined with O displacementsQ'U, and on-sit[e-|atp 
tractive interaction U based on non-linear screeningt2lllil, 
have been proposed as the origin of superconductivity in 
doped BaBiOa. In these approaches, repulsive interac- 
tions at finite distances play an important role, what is 
consistent with poor screening of Coulomb interactions 
by the low density of carriers. |— . 

As first shown by Rice and SneddontS, treating the dis- 
placements d of O ions in the direction of their nearest- 
neighbor Bi ions in the antiadiabatic approximation, 
leads to a decrease in U by zg^ /K, and an increase of the 
repulsion V between nearest Bi atoms by /K, where z 
is the coordination number, g the electron-phonon inter- 
action, and K the second derivative of the elastic energy 
with respect to d. This leads naturally to the extended 
Hubbard model for the description of doped BaBiOa, and 
to superconductivity carried by bipolarons (a pair oLcac.- 
riers at Bi sites accompanied by O displacements)E30. 
While the excitonic mechanism involves partial occu- 
pation of O statfisJpy holesU, in agreement with opti- 
cal experimentsEililj, it is still possible that this one- 
band modisl describes the low-energy physics (as in the 
cupratesEZl) . This is certainly the case in the purely elec- 



tronic model proposed by Varma, in which non-linear 
screening not only reduces the bare atomic U ^ 11 eV, 
but renders it negatival3. 

In recent years, the phase diagram of the extended 
Hubbard model in one dimension has been studied by nu- 
merical techniques, with particular emphasis in the|«HiM|- 
ter filled case (number of particles per site n = l/2)ll3"c3. 
For n = 1/2 or n — 2/3, V > 5 and small negative 
U, calculations of the correlation exponent Kp in sys- 
tems of up to L = 16 sites, predict a Luttinger liq- 
uid phase with dcimijiant superconducting correlations 
at large distancealj'Ej. However, from an analysis based 
on the infinite V limit and the extrapolation of the spin 
gap, Penc and Mila suggested that the results of Kp were 
affected by finite-size effects, and the_aystem would not 
be a Luttinger liquid in that regionlij. For n = 1/2, 
perturbative arguments around the infinite V limit (de- 
veloped in more detail here) suggested that this region 
corresponds to phase separationcj. Monte Carlo studies 
in systems with L ^ 64 lattice sitesE^ have shown that 
the region of phase separation in the model extends to 
much lower values of V than those obtained for L < 16, 
leaving practically no place for a Luttinger liquid phase 
with dominant superconducting correlations for large V. 
The extent of this phase remains unclear. Furthermore, 
due to the above mentioned finite-size effects, and tejchni- 
cal problems with different Monte Carlo methodsS, the 
region of small \U\ and V > 8t, where t is the hopping, 
is practically unaccessible to the present available nu- 
merical methods. This region is also out of the range 
of applicability of the continuum-limit field theory (also 
called g-ology), which has beeripapplicd to the extended 
Hubbard and related modelsc3ca. 

In this work, we study the one-dimensional attractive 
Hubha^rd model in the limit V >> t. Extending previous 
workll3, we derive an effective low-energy Hamiltonian 
including terms up to second order in t. Under certain 
approximations, which are not essential for sufficiently 
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large V, the effective Hamiltonian is mapped into two 
XXZ models, representing the movements of particles 
in singly and doubly occupied sites respectively. From 
the Bethe ansatz solution of these models, the phase di- 
agram in the thermodynamic limit L — *■ cxd as a function 
of n, U/t < and large V/t is obtained. In the next Sec- 
tion, we explain the model and the mapping procedures. 
The results are shown in Section III. Section IV contains 
the conclusions and a discussion of the possible relation 
between our results and the phase diagram measured in 
Bai.^K.BiOaQ. 

II. MODEL AND EFFECTIVE HAMILTONIAN 

The one-dimensional extended Hubbard Hamiltonian 
in standard notation is: 

H = -t^{4+icyCia + H.c.) + U^Ui^riii + rii+irii, 

i(7 i i 

(1) 

and we consider the case U < V » t. In this 
limit, there are no nearest-neighbor occupied sites in 
the low-energy subspace. The terms in Eq. which 
mix states of this subspace with states with non-zero 
occupancy at nearest-neighbor sites, can be eliminated 
through a staiadard canonical transformation, as de- 
scribed in Ref.tJ. This procedure originates terms of 
order /V or smaller in the low-energy subspace. The 
resulting effective Hamiltonian within this subspace (a 
projector P = Wi^a'i^ ~ n-ia-rii^icr') is implicit) can be 
written in the form: 

^eff — + Hd + Hsd, (2) 
where Hg involves only singly occupied sites: 

Hs = -t^icl+la^ta +H.C.)(1 -n^+lff) 
ia 

+ Vs ^ nicrni+2a' (1 - ni^){l - ni+2a) 

+ VsJ2'^4+3a('i+2a4+la'''-^a' + H.C.), (3) 
iaa' 

with Vs = —t'^/V. The second term is an interaction 
between second nearest neighbors, and the third term 
displaces two next-nearest-neighbor particles one lattice 
parameter. 

The second term in Eq. (^) can be described entirely 
in terms of the operators ~ c^jC^j^j which create "dou- 
blons" at each site: 

Hd = AdJ2 4d^ - td Y^idl+A + H.c.) 

i i<7 

+VdY,4^^d,+2dld,, (4) 



where Ad = U - At^/iV - U), td = 2tV(V - U), and 
Vd = 4£^[1/{V -U)- 1/{W - U)] are the effective on- 
site energy, nearest-neighbor hopping, and next-nearest- 
neighbor repulsion respectively, for doublons. 

The last term in Eq. (^) describes interactions between 
singly and doubly occupied sites: 

Hsd = Ml Y.[{c\^^^c\_,^ + cLi^ct^i^) 

i 

x(l - nj_i)(l - nj+i)(di_i + 2di + di+i) + H.c] 

+M2 Y,{4ad\+2d^C^+2a + H.C.) 

+Vsd ^ nia{l-nis)d\j^^di+s. (5) 

ia&=±2 

Here Mi = -t^[l/V + 1/{V - U)]/2 describes annihila- 
tion of a doublon with creation of two particles at empty 
sites, and the Hermitian conjugate process, while M2 = 
e/{2V - U) {Vsd = t'[2/(T/ -U)- 2/{2V- U) - \/{2V]) 
corresponds to interchange (interaction) of a doublon and 
a particle at a singly occupied next-nearest-neighbor site. 

For V = -fcx), all terms of i?e// vanish except the first 
term of Eq. (^) and the first term of Eq. (jj) , and H^f / 
can be solved exactlyllj. In this limit, for U > ~4t and 
sufficiently small density n, the system has no doubly 
occupied sites and is described by Hs- Instead, for U < 
—At, n < 1, the ground state of H^ff is the same as 
that of Hd and has no singly occupied sites. For other 
values of U and n < 1, the system phase separates into 
the phases just described, and the limits of the region of 
phase separation (PS) can be obtained using the Maxwell 
construction (finding the common tangent to the curves 
Es{n) and Ed{n), where Ea{n) is the ground state of 
at density n). 

For the sake of clarity we call "metallic" (M) the phase 
without double occupancy (ground state of Hs), and 
"bipolaronic" (BP), the phase described by the ground 
state of Hd, although the negative U is not necessarily 
related with atomic displacements in a real system. For 
finite but large V, the energy cost for constructing a uni- 
form phase with both, singly and doubly occupied sites 
is high in comparison with Mi, M2 and Vsd- Thus, in the 
ground state, Hsd can only act in the PS region, at the 
boundary between M and BP phases, and is irrelevant 
in the thermodynamic limit. Our main approximation 
is the neglect of Hsd- This should be correct as long as 
the energy gain of Hs in the PS region (~ t) is larger 
than the terms of Hsd (~ ^^/^)- We also neglect the 
last term of Eq. (^) . This term vanishes at the extreme 
densities of the M phase (n = and n = 1/2). The ratio 
of its expectation value with respect to the expectation 
value of the first term of Eq. (|^) can be estimated by 
perturbation theory: 

r = -lr[nsCOs{Trns) - -sm{Trns)], (6) 
V n 

where — n/(l — n). The remaining terms of Hs, 
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and the whole of Hd can both be mapped into a spin- 
less fermion model H^^ {a = s or d): 

Hit = A„ ^ - Y.^fl,h + H.c.) 

+Vjlfjlj,+u (7) 

with As = 0. When a = s, a site occupied by a fermion 
fl corresponds to a single occupied site and an empty 
site aLthe right of it, as explained in detail by Penc and 
MilaE£l (a similar mapping was also-used in a model for 
oxygen ordering in YBa2Cu306+a:L3). Then, for L sites 
and N particles in Hg , the corresponding number of sites 
and particles in H^^ are = L — N, = N. Then, 
the energy per site es{n) of Hg for density n is related to 
the corresponding quantity e''/{ns) of H^-^ by: 



(8) 



Similarly, can be cast into the form of Eq. , map- 
ping a doubly occupied site and an empty site at the 
right of it into a single site occupied by a fermion. The 
mapping of the different physical quantities is the same 
as that used before to find the correlation exponent Kp 
in a generalized t — J model with very large three-site 
tcrmEJ'Eil. The number of sites and fermions in H'^^ are 
Ld = L- N/2, Nd = N/2. The energy per site of Hd is 
given in terms of that of H^J by: 



ed{n) = (1 - ^)ej{nd)\ 



nd 



(9) 



To calculate Kp, we also need the mapping of the 
velocityL2l: 



L sf 



and/or the Drude weight: 

Ld^Ed^ _Ldd^Elf{<^sf) 



L 



(10) 

2(2 -n)Z?y. 

(11) 



Here Ed{<i>) {E'/{<^sf)) is the energy of a ring described 
by Hd (Hf) threaded by a flu^p («>«/)■ The correlation 
exponent can be calculated astil: 



TTVd 



2d'^ed/dn? Vd 



Dd 



2d'^ed/dn- 



1/2 



(12) 



Similar expressions give Kp for Hs, but we do not give 
them, since in the M phase always Kp < 1, and we are in- 
terested in the region Kp > 1, for which superconducting 
correlations dominate at large distances. 

Using a Jordan- Wigner transformation, H^^ is trans- 
formed into an equivalent XX Z model with La sites and 
Ma spins down: 



HfJ = ijf^^ = 2U 



{Sf S: 



qy qV 



Va 



+ AaNa + VaiNa~ La/i). 



(13) 



We have calculated the energy e'^^ [ua) solving numer- 
ically the integral equations of the exact Bethe jansatz 
solution of Eq. (^3|) in the thermodynamic limito. To 
obtain Kp, we have calculated the excitation energies for 
two small momenta, solvingjijimerically the correspond- 
ing Bethe ansatz equation^. This allowed us to ex- 
tract v''^ {ud). From it, the numerical second derivative 
of ed(n), Eqs. (|), ^ and the first Eq. (|2|), Kp was 
calculated for Ud ^ 1/2. For Ud = 1/2, we had technical 
problems in the calculation of v'^J-Lori), but fortunately, 
analytical expressions are knownOEJ: 



f _ TTtsin^ f 

"d — ,, ' ^d 



^ = a.Tccos{V/2t). 



(14) 



Using Eqs. (|ol), ([III), (|ll|) and (|T|), one has for n = 2/3 
(ud = 1/2): 



Kpi2/3) 



Hl-fi/n)- 



(15) 



For Va = 0, H^f can be solved triviallyEl and Eqs.(|), 
( |lO| ) , ([l^) , ( [l^ ) lead to another analytical result for Vd = 



ifp(2/3) 



(2 - nf 



(16) 



Since always Vrf/i^ > 0, and Kp decreases with increasing 
Vd, Eq. ( [iq ) implies that for large V, no phase with 
dominant superconducting correlations exists for n> 2 — 
\/2 ~ 0.59. 



III. BETHE ANSATZ RESULTS 




FIG. 1. Ground state energy of the phase with singly oc- 
cupied sites (dashed line) and doubly occupied sites (full line) 
as a function of density. The dotted line is the Maxwell con- 
struction (see text). 
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In Fig. 1 we show the energy per site of Hg (es(n)) 
with the last term neglected, and that of Hd (ed(n)), as 
a function of density n for [/ = —2 and V — 7. We take 
i = 1 as the unit of energy. The dotted line represents 
the energy for average composition n, of an inhomoge- 
neous (phase separated) mixture of the "metallic" (M) 
phase described by Hs for density ni — 0.2307 and the 
"bipolaron" (BP) phase (ground state of Hd) for den- 
sity ri2 = 0.5901. These compositions nt, represented 
by diamonds in Fig. 1, are obtained finding the com- 
mon tangent to both curves es(n) and ed{n) (Maxwell 
construction). Between them, the energy of the phase 
separated phase is lower than both es(n) and ed{n). The 
energy es{n) is dominated by the first term of Eq. (||), 
which already exists for V = +oo. The effect of finite 
large V is small (except near n — 1/2 for which the first 
term of Eq. (|^) vanishes), and does not change ni sig- 
nificantly. Instead, the effect of a finite large V on n2 is 
dramatic, reducing it from 712 = 1 to n2 < 0.6. This is 
because for V ~ +00, ed{n) is always a straight line (ex- 
tending from ed(0) = to 6^(1) = -1 ii U = -2). The 
second and third term of Hd (Eq.(^)), taken into account 
exactly, are responsible of the curvature of ed{n) and the 
shift of 712 from 1 . This is the main effect of finite large 
V. 



K1.5 - 




FIG. 2. Correlation exponent Kp as a function of density 
in the "bipolaron" (BP) phase for t/ = -2 and 1/ = 7 (full 
line). The dashed line is the result in absence of the effec- 
tive nearest-neighbor repulsion between doubly occupied sites 
(V, = 0). 



on U and V within the studied range of parameters, de- 
creasing with increasing repulsions. For U = —2, V = 7 
(as in Figs. 1 and 2), ric = 0.429. These results disagree 
with tho s ii o btained by numerical diagonalization of small 
systemsE3E3, which obtained Kp > 1 for small values of 
\U\, large V and n — 1/2, n — 2/3, but agree with the 
statement that the^ system might not be a Luttinger liq- 
uid in that regiontj, and with recent Monte Carlo results 
in larger systemsa. These results and the ones shown be- 
low indicate that there is phase separation (PS) in that 
region. We believe that the reason of the artificially large 
values of Kp in the above mentioned numerical results is 
that they were calculated using the first Eq. ( p^ ) with 
d^e/dn? ~ l/(Kri^) determined numerically from the en- 
ergy for N, N — 2 and TV -I- 2 particles with N/L = n. 
In phase separated regions (dashed line in Fig. 1), the 
compressibility k diverges in the thermodynamic limit, 
but in small systems, k can be large and positive due to 
finite-size effects, leading to very large values of Kp, while 
in fact the system is not a Luttinger liquid. This effect 
was present in numerical studies of tlie one-dimensional 
t — J model with correlated hoppingpj. 
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Another important effect of a finite V is that the pairs 
acquire mobility and for low densities, superconducting 
correlations dominate at large distances [Kp > 1) in the 
BP phase. In Fig. 2 we represent Kp(n) in this phase for 
the same parameters of Fig. 1. Also shown is the analyt- 
ical result Eq. (p^, which is an upper bound of Kp for 
any values of C/ < and V » t. For V -l-oo, we ob- 
tain a critical density Uc — 0.422 for which Kp(nc) — 1. 
For n < Tic and arbitrary values of ?7 < and V >> t, 
the exponent Kp[n) > 1. Thus, ric lies in the interval 
(0.422,2 - V2). For -5 < [/ < and 5 < F < 20, 
we find that Uc < 0.45. Also, ric depends very weakly 



FIG. 3. Phase diagram of the model in the density-f/ plane, 
for different values of V . The phase without doubly occupied 
sites is called "metallic" (M), while the phase with no singly 
occupied sites is denoted as "bipolaronic" (BP), and if Kp > 1 
we call it "superconducting" (S). The region of phase separa- 
tion is labeled PS. 
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In Fig. 3 we show the phase diagram of the model, 
determined using the MaxweU construction (as in Fig. 
1) for different values of U and V (points indicated with 
solid symbols in Fig. 3), and searching the critical den- 
sities He for which Kp{nc) = 1. The phase diagram for 

V — > +00 was already known, except for the boundary 
a,t n — 0.422, which separates the BP phase with Kp < 1 
(in which charge correlation functions are the dominant 
ones at large distances) from the "superconducting" (S) 
phase with Kp > 1 . The most noticeable effect of a finite 
V, already present for V as large as 20t is the change 
in the boundary between the BP and PS regions, as a 
consequence of the above mentioned curvature of ed{n) 
for V < +00. Also, as V decreases, the region S with 
Kp > 1 increases, moving to larger values of U and to 
slightly larger densities. Finally, the PS region is reduced. 

The case V = 5is probably beyond the quantitative va- 
lidity of our large V approximation. Monte Carlo results 
for V = 8 (Fig. 14 of Ref.Ej) are in qualitative agreement 
with our results. Quantitatively, Clay et al. obtain that 
for U = 0, phase separation begins at ni ~ 0.5, while 
ior V — 7 we obtain ni = 0.345. This difference might 
be due to the effect of terms of order t"^ with m > 2 or 
the third term of Hg (Eq. (^)), which we have neglected. 
Terms of order t'* reduce td and increase ni. The limit 
between S and M phases at low densities is affected nei- 
ther by Hsd (Eq. (||)) nor by the third term of Hs, but 
might be changed slightly by terms of order i^. 

IV. SUMMARY AND DISCUSSION 

Generalizing a previous approachEl, we have con- 
structed an effective Hamiltonian i?e// for the extended 
Hubbard model when V >> t and U < 0. This permits 
to study the region V > 8t, which is outsidc-|thc re- 
gion of validity of weak-coupling approaches^ .ca^while 
numerical diagonalization of small systemaiaEJ'ES dis- 
play important finite-size effects, as discussed in Section 
III, and Monte Carlo calculations have some technical 
problemso. The effective Hamiltonian can be divided in 
three parts: Hs, Hd and Hgd- The latter is irrelevant for 
sufficiently large V , and Hs {Hd) acts on a phase in which 
all particles move in singly (doubly) occupied sites. Hd 
and the most important terms of Hs were mapped into a 
Bethe ansatz exactly solvable model, which allows to ob- 
tain the energy and correlation exponent Kp in the ther- 
modynamic limit. From this information we have con- 
structed the phase diagram. We obtain a region for low 
and intermediate densities and sufficiently negative U , in 
which the system behaves as a Luttinger liquid with dom- 
inant superconducting correlations at large distances, for 
any finite V » t. The existence of this phase and the 
main changes in the phase diagram with respect to the 

V = +O0 limit, are due to the dynamics of doubly occu- 
pied sites in Hd, controled by terms of order /{V ~U). 
The most noticeable effect of a finite V for moderate val- 



ues of \U\ {U > —4), is that the upper density of the 
phase separation (PS) region is reduced from n = 1 to 
much lower values. Instead, the larger density for which 
superconducting correlations dominate at large distances 
remains near n — 0.42. Thus, the place left by the PS 
region is mostly occupied by the BP ( "bipolaronic" "nor- 
mal") phase rather than by the S ("superconducting ") 
phase. Nevertheless, as V decreases from very large val- 
ues toV^7, the upper value of U for which the S phase 
exists increases from -4 to ~ —3. 

The extended Hubbard model with attractive U can be 
justified in different ways, as a model for doped BaBiOa, 
as discussed in Section I. In spite of the different dimen- 
sionality of the real compound, one can discuss qualita- 
tively the physics expected from the phase diagram (Fig. 
3). BaBiOs has one electron per site (n = 1), and the 
ground state of the model is a charge density wave, in 
which nearest-neighbor .sites are not equivalent, as ex- 
perimentally observedQT]. For sufficiently negative U, as 
n decreases (corresponding to partial replacement of Ba 
for K) keeping ric ~ 0.4 < n < 1, the ground state of 
the model is a Luttinger liquid with dominant charge 
density wave correlations at large distances, with wave 
vector 2kp = nn. Experimentally, the ground state of 
Bai_2,Ka;Bi03 has different charge density wave order- 
ings for < a; <'-^ 0.4 (1 > n >~ 0.6)0. As n is fur- 
ther lowered (corresponding to increasing x = 1 — n), 
the model enters a region with dominant superconduct- 
ing correlations at large distances, which has a corre- 
sponding superconducting phase in Bai-ajK^jBiOa (for 
^ 0.4 < a; < 0.5)lI. Above x — 0.5 the experimental 
system cannot be formed, since the solubility limit of K 
atoms is exceeded. This picture is also consistent with 
mean- field calculations rin the three-dimensional model 
for different parameterai3. 
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